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ABSTRACT 



>< 



Aims. Previous observations with the H.E.S.S. telescope array revealed the existence of extended very-high-energy (VHE; 
1 E>100 GeV) 7-ray emission, HESS J1023-575, coincident with the young stellar cluster Westerlund 2. At the time of discov- 
ery, the origin of the observed emission was not unambiguously identified, and follow-up observations have been performed to 
further investigate the nature of this 7-ray source. 

Methods. The Carina region towards the open cluster Westerlund 2 has been re-observed, increasing the total exposure to 45.9 h. 
The combined dataset includes 33 h of new data and now permits a search for energy-dependent morphology and detailed spec- 
troscopy. 

Results. A new, hard spectrum VHE 7-ray source, HESS J1026-582, was discovered with a statistical significance of 7a. It is 
positionally coincident with the Fermi LAT pulsar PSR J1028-5819. The positional coincidence and radio/7-ray characteristics of 
the LAT pulsar favors a scenario where the TeV emission originates from a pulsar wind nebula. The nature of HESS J1023-575 
is discussed in light of the deep H.E.S.S. observations and recent multi-wavelength discoveries, including the Fermi LAT pulsar 
PSR J1022-5746 and giant molecular clouds in the region. Despite the improved VHE dataset, a clear identification of the object 
responsible for the VHE emission from HESS J1023-575 is not yet possible, and contribution from the nearby high-energy pulsar 
and/or the open cluster remains a possibility. 



Key words. Gamma-ray: stars; ISM: HII regions - ISM: individual objects: HESS J1023-575, HESS J1026-582, RCW 49 (NGC 
3247, G284.3-0.3), PSR J1028-5819, PSR J 1022-5746 



1. Introduction 

The H.E.S.S. (High Energy Stereoscopy System) collabora- 
tion reported the detection of a bright, extended (standard 
deviation g=0.18°) VH E 7-ray source HESS J1023-575 
(|Aharonian et al]|2007a|) coincident with the open stellar 
cluster Wcstcrlund 2. This discovery provides evidence 
for particle acceleration to extremely high energies in 
this region. Young stellar clusters are among the sites 
considered for high-energy 7-ray production. This relates 
to the existence of late states of stellar evolution, which 
are generally considered favorable for particle acceleration. 
Different scenarios have been proposed to account for 
VHE 7-ray emission from stellar environments, consid- 
ering the contribution of hot and massive Wolf-Rayet 
(WR) and OB stars and their winds in single, binary, 
or collective processes, pulsars and their synchrotron 
nebulae, as well as supernova explosion s and their ex- 
pand in g remnants (Eichler fc Usovl 119931: IWhite fe Chenl 
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19951: iBenaglia fc Romero! 120031: iReimer et al.l 120061 



Pittard fe Dougherty! 12001 iTorres et al.l 12004 iBednarekl 



20051: iManolakou et all l2007t) . The giant HII region 
RCW 49 (NGC 3247) and its ionizing cluster Westerlund 2, 
hosting over two dozen of massive stars as well as two 
remarkable WR stars (WR 20a and WR 20b), have 
been studied extensively over the whole electromagnetic 
spectrum, and certainly interest was further boosted by 
the discovery of VHE 7-ray emission from the vicinity of 
Westerlund 2. 

The VHE 7-ray source HESS J1023-575 
(jAharonian et al.l l2007aT) is characterized by a power- 
law spectrum with a photon index T w 2.5 and, at a 
nominal distance of 8 kpc (ob tained from timing studies 
of the WR binary WR 20a; iRauw et al.l I2007T) . would 
have a 7-ray luminosity L 7 (E > 380 GeV) = 1.5 xlO 35 
erg s _1 . The extent of the observed VHE 7-ray emission 
spatially overlaps with the wi nd-blown bubbles obs erved 
in infrared (IR) by Spitzer (jWhitnev et all l2004f ) and 
in radio continuum by A TCA at 1.38 and 2.38 GHz 
(jWhiteoak fc Uchidal 1 19971) . One of these bubbles, with 
a diameter of —7.3' is believed to be connected to the 
massive binary WR 20a, surrounding the central region 
of Westerlund 2, while the second has been related to 
WR 20b, with a diameter of -4.1'. The WR binaries 
WR 20a and WR 20b have been dete cted in X-rays in 
the 0.5 to 10 keV band with Chandra (jNaze et all I2008D 
with luminosities of 19xl0 32 erg s -1 and 52xl0 32 erg 
s _1 respectively. The high X-ray luminosity seen from 
WR 20a, which is likely associated with the open cluster, 
provides support for the scenario whereby in a colliding 
wind zone energ y dissipated in s hocks is converted into 
thermal energy (|Naze et all 120081 ). The observed source 
extension (0.18°) and non-variability of the 7-ray flux 
of HESS J1023-575 disfavor any scenario where particle 
acceleration occurs close to the massive stars, and in the 
colliding wind zone within the massive WR 20a binary 
system scenario. A pro minent feature at radio wave lengths, 
known as the "blister" (IWhiteoak fc Uchidalll997lh coinci- 
dent with HESS J1023-575 (jAharonian et al.ll2007al) . does 
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not seem to be evident in the diffuse X-ray flux measured 
by Chandra (|Naze et al.ll2008l ). which appears featureless 
and predominantly of thermal origin. Upper limits on the 
non-thermal diffus e emission were recently reported by 
iFuiita et al.l (|2009l ) in the 0.7-2 keV (F ul <2.6xl0- 12 erg 
cm" 2 s" 1 ) and 0.7 to 10 keV (F ul <6.1 x 10~ 12 erg enr 2 
s _1 ) bands, using Suzaku observations and integrating 
ove r the entir e came ra field (17.8' x 17.8'). According 
to IFuiita et al.l ((2009) , the statistics are not sufficient 
to distinguish between three thermal components and 
two thermal components plus a non-thermal one when 
modeling the background. They claim that their mul- 
ticomponent fit supports the existence of non-thermal 
emission. However, the ambiguity in the goodness-of-fit 
using a multi-thermal component and a multi-thermal 
component including a non-thermal one does not warrant 
such a conclusion on the basis of the presented data. 

A possible scenario discussed to explain the VHE 7-ray 
emission includes the interaction of cosmic-ray protons ac- 
celerated in stellar wind shoc ks or supernova blast wave s 
with dense molecular clouds ([Aharonian fc Atovanlll996h . 
The molecular cloud content in this region observed by the 
CfA 1.2 m Millimeter- Wave Telescope and the 4-m Nantcn 
2 sub-mm tel escope shows a complex and rich environment. 
iDamd (|2007l) observed Westerlund 2 with the CfA 1.2 m 
telescope using 12 CO (J = 1-0) survey data and found a 
possible association with a giant molecular cloud (GMC) 
on the far side of the Carina spiral arm, at a distance of 5.0 
kpc. Observations at a higher resolution (90") on the J = 2- 
1 transition of 12 CO performed by the Nanten 2 4-m sub- 
millim eter telescop e dFurukawa et al.l 120091 : lOhama et al.l 
I2010D confirm the IDamd (|2007t) results but also suggest 
another possible association. A significant clump of molec- 
ular gas, consisting of a foreground component at km s~ 1 
is found, with a mass of (9.1±4.1) x 10 4 M , and a back- 
ground component at 16 km s _1 , with a mass of (8.1±3.7) 
xlO 4 Mq, both of which are claimed to be connected with 
RCW 49, at a distance of 5.4ti;4 kpc. The analysis of the 
Nantcn survey in the 12 CO (J = 1-0) (|Fukui et al.ll2009l) 
also shows two unusual structures in the vicinity of the line- 
of-sight projection: a straight feature in the East and an 
arc-like feature in the West in the velocity range 20-30 km 
s _1 , corresp onding to a distance of —7 kpc (using the rota- 
tion curve of lBrand fc Blitz|ll993l) , although the association 
with the open cluster is unclear. 

At high energies (HE; E>100 MeV) the Fermi Large 
Area Telescope (LAT) collaboration has recently reported 
the detection of two new sources: 1FGL J1023. 0-5746, posi- 
tionally coincident with HESS J1023-575, and a new 7-ray 
pulsar 1FGL J1028.4-5810 (PSR J1028-5819) to the east of 
the H.E.S.S. sour ce and associated w ith the EGRET source 
3EG J1027-5817 (lAbdo et al.ll2009al) . 

The first source, reported in the bright source list, 
has been very rec ently identified with a new 7-ray pulsar , 
PSRJ1022-5746 (iDormodv et al.l 120091 lAbdo et all 120101: 
ISaz Parkinson et ah l2010h . The pulsar periodicity found 
in the blind search performed by the Fermi LAT collab- 
oration is 111.47 ms. The young pulsar (with a charac- 
teristic age of r —4.6 kyr) has a spin-down luminosity of 
1 . 1 x 10 37 erg s~ 1 . Radio pulsations have not been found but 
the 130 ksec Chand ra image shows a faint source (CXOU 
J102302.8-574606.9. lTsuiimoto et al~ll2007l ) identified as a 
possible counterpart. The second source, PSR J1028-5819 



H.E.S.S. Collaboration: Revisiting the Westerlund 2 Field with the H.E.S.S. Telescope Array 



3 




(|Abdo et al.ll2009al Manchester et alJl2005h is a 91.4 ms 7- 
ray and radio pulsar with a total spin-down power of E 
= 8.3 xlO 35 erg s _1 and a phase-average luminosity from 
0.1 to 30 GeV of (120±73)xl0 33 erg s" 1 at a dispersion -57°00 
measure-derived distance of 2.3±0.7 kpc. Its characteris- 
tic age is 8.9 xlO 4 yrs. The HE photon spectrum shows a 
hard cutoff at 1.9±0.5 GeV, while no contribution of steady 
emission is reported. s 
The unresolved origin of the VHE 7-ray emission from 8 
the direction of Westerlund 2 motivated further observa- ^ 
tions, aiming to discriminate among the alternative scenar- a 
ios for the VHE emission. The results of these deep obser- ^ -58°00 
vations, alongside the discovery of a new, hard, VHE 7-ray a 
emission coincident with the Fermi LAT source PSR J1028- 
5819, are reported here. The nature of the VHE sources is 
then discussed in light of the new multi-wavelength obser- 
vations in the Westerlund 2 field. 



2. Analysis and Results 

The High Energy Stereoscopic System (H.E.S.S.) is an ar- 
ray of four VHE 7-ray imaging atmospheric Chercnkov 
telescopes (IACTs) located in the Khomas Highland of 
Namibia (23°16'18"S 16°30'00" E). Each of these telescopes 
is equipped with a tessellated spherical mirror of 107 m 
area and a camera comprised of 960 photomultiplier tubes, 
covering a large field-of-view (FoV) of 5° diameter. Its large 
FoV and good off-axis sensitivity make H.E.S.S. ideally 
suited for studies of extended sources and regions of the sky 
in w hich more than one 7-r ay source could be present (see 
e. g.[Ah aroman et al.ll2008D The system works m a comci- 
dence mode (e.g. lFunk et al.ll2004h . requiring at least two of 
the four telescopes to trigger the detection of an extended 
air shower (EAS). This stereoscopic approach results in a 
high angular resolution of ~ 5' per event, good energy res- 
olution (6% on average) and a n effective background re- 
jection (jAharonian et al.ll2006bh . These characteristics al- 
low H.E.S.S. to reach a sensitivity ~2.0xl0~ 13 ph cm _2 s _1 
(equivalent to 1% of the Crab Nebula flux above 1 TeV), 
or less if advanced tech niques are used for image analysis 
(|Aharonian et al.ll2009d) , for a point-like source detected at 
a significance of 5cr in 25 hours of observation at zenith. 

Previous observations of the re gion containing 
PSRJ1028-5819 and Westerlund 2 (jAharonian et al.l 
2007a) were performed between March and July 2006, 
amounting to a total of 12.9 h of live time. The analysis of 
the data led to the discovery of an extended source (a = 
0.18°±0.02°) located at a=10 h 23 m 18 s and <5=-57°45'50" 
(J2000). The morphological and spectral analyses of the 
source did not allow a firm identification of the origin 
of the VHE emission. Energy-dependent morphology 
studies have proven to be a powerful tool to distinguish 
between alternative source models and respective VHE 
7-ray emission signatures (|Aharonian et al.1 l2006al) . In 
order to discriminate among the different 7-ray emission 
scenarios proposed and to investigate energy-dependent 
morphology, follow-up observations were pursued. The 
region was observed for an additional 33 h in January 
and March 2007, April and May 2008, and May to June 
2009. The complete data set is presented here and has 
a live time of 45.9 h after quality cuts to remove data 
taken du ring bad weather condit ions or hardware irreg- 
ularities (jAharonian et al.l l2006bl ). The data were taken 
in wobble-mode configuration, where the telescopes are 
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Fig. 1: Image of the Westerlund 2 field showing the significance 
of the VHE 7-ray emission using an oversampling radius of 0.22°. 
The contours in black (dashed) and green (solid) show the sig- 
nificance levels for the maps at low and high energy range, re- 
spectively, using the same oversampling radius, starting at 5a 
in steps of la. The maximum of the color scale has been chosen 
such that the typical detection threshold ~5cr appears on the 
transition from blue to red. 



pointed offset (from 0.5° to 2.5° in this data set) from the 
nominal source location to allow simultaneous background 
estimation. The zenith angles range from 30° to 55°, 
with a mean value of 38°, which results in a mean energy 
threshold of ^0.8 TeV when using hard cuts on the EAS 
image size of 200 photoelectrons (p.e.). The mean energy 
threshold of the new data set is slightly higher than the 
one in the previous analysis due to lower reflectivity of the 
IACT mirrors compared to the 2006 observations and due 
to a stricter cut on the image size than the one used for 
the 2006 data set (standard cuts at 80 p.e.). The optical 
response of the system was estimated from the Cherenkov 
light of single muons hit ting the telescopes as explained in 
lAharonian et al.l (|2006blh 

The data have been analyzed using the H.E.S.S. stan- 
dard event reconstruction scheme (jAharonian et ah] 2006bfl 
using the Hillas second moment method (iHillasI 119850 
second independent analysis chain was used to cross-check 
the results, including independent calibration of pixel am- 
plitudes and identification of problematic or dead pixels 
in the IACTs cameras, leading to compatible results. To 
achieve a maximum signal-to-noise ratio and to improve 
angular resolution, hard cuts were used to reduce the data, 
and only observations runs (of 28 min duration) with four 
active IACTs were used to study the morphology and spec- 
trum of the source. 

The large energy range covered by H.E.S.S. allows the 
investigation of the morphology and spectral features in 
the region of interest in different energy bands. For this 
analysis, two energy bands were selected a priori such that 
each band has a similar signal-to-noise ratio, for the range 
of zenith angle of the current observations, assuming the 
source has a photon index of 2.2 (a typical value for Galactic 
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Fig. 2: Gaussian-smoothed (cr=0.08°) excess images for the Westerlund 2/RCW 49 and PSR J1028-5819 region. On the left, the low 
energy map (0.7 TeV <E < 2.5 TeV) is shown, while the high energy map (E > 2.5 TeV ) is displayed on the right. The H.E.S.S. 
significance contours (dashed black lines for the low energy map and solid green ones for the high energy map) are calculated using 
an oversampling radius of 0.1° and are shown above 4a in steps of la. The position of WR 20a and the two Fermi LAT pulsars 
are marked in cyan and white (for more details see text). 



VHE 7-ray sources). Images of the VHE 7-ray excess (see 
next Sect.), corrected by the corresponding radial accep- 
tance in each energy band, were produced for two bands 
(0.7< E<2.5 TeV and E>2.5 TeV). 

2.1. Energy-Dependent Morphology 

Fig. Q] shows the significance map of the region using 
an over sampling radius of 0.22 ° (optimized for extended 
sources. [Aharonian et al.l l2006ch . The background in each 
bin of the image was es timated using a ring with a mean ra- 
dius of 1° (as defined in lAharonian et al~ll2006bl) around the 
test position. The previously reported source HESS J1023- 
575 is detected with a peak significance of 16a corre- 
sponding to 545 excess events above background. A sec- 
ond excess appears to the south-east towards the direction 
of PSRJ1028-5819 (dubbed HESS J1026-582), showing a 
peak significance of 7a pre-trials (169 excess events). Source 
locations are extracted from the uncorrelated excess map 
fitted to a two-dimensional Gaussian function folded with 
the H.E.S.S. point spread function (PSF; -0.08° above 200 
p.e.). The best fit position lies at a = 10 h 23 m 24 s ±7.2| tat 
and 6 = -57°47'24"±l'12" stat (J2000) for HESS J1023- 
575, compatible with the previously reported position. For 
the second source, HESS J1026-582, the Gaussian peaks at 
a = 10 h 26 m 38.4 s ±21.6| tat and S = -58 o 12'00"±l'48" st at 
(J2000), 0.6° apart from the other so urce. The syste matic 
error is estimated to be 20" per axis (|Gillessenll200l . The 
two derived positions and statistical errors are marked with 
black and green crosses in Fig. [TJ 

To further investigate the nature of the two detected 
excesses, the image has been analyzed in the two energy 
bands described above. For the lower energy band, between 
0.7 and 2.5 TeV, the image (Fig. [5]) shows a distinct ex- 
cess emission at the location of HESS J1023-575, extended 
beyond the nominal size of the H.E.S.S. instrument PSF, 
with an intrinsic extension of er=0.18 o ±0.02°. The second 
excess is strongly reduced in this energy band. However, 



this second emission region, HESS J1026-582, is clearly vis- 
ible on the high energy image (E>2.5 TeV), indicative of 
a very hard photon index and also evidence for extension 
at the scale ct=0.14°±0.03°. The two Gaussian-smoothed 
(cr=0.08 o ) excess maps are shown in Fig. [2] with signifi- 
cance contours above 4er in black dashed lines for the low 
energy map and solid green lines for the high energy one, 
obtained by using an integration radius of 0.1°. The signifi- 
cance contours for a larger integration radius (0.22°) for the 
lower (in black) and the higher (in green) energy bands and 
are also shown in Fig. [1] together with the total significance 
map above 200 p.e. 

To further investigate the interplay between the 
two neighboring sources as of HESS J1023-575 and 
HESS J1026-582, a rectangular region (slice) is defined in 
the uncorrelated excess maps along the line connecting the 
best-fitted positions of the two spots with a width of twice 
the H.E.S.S. PSF (see Fig. |3] inset). The profile resulting 
from a projection of the slice on the long axis is shown 
in Fig. [3] for the low (in red) and high (in blue) energy 
images. The radial profile for the low energy map is well 
fitted by a single Gaussian function (x 2 /^=5.57/6, corre- 
sponding to a probability of P=0.47) centered at a position 
compatible with the centroid of HESS J1023-575 and the 
position of Westerlund 2 and PSR J1022-5746 (cyan and 
white markers respectively in the inset figure). Their po- 
sitions with respect to the center of the slice are marked 
with dashed lines in the profile figure. However, the pro- 
file corresponding to the high energy events shows a sec- 
ond peak towards the direction of PSR J1028-5819 (dashed 
line) and a fit to a double Gaussian function is clearly fa- 
vored (x 2 /^=0.45/4, P=0.97) against a single Gaussian one 
( x 2 /^=10.88/6, P=0.09). 

Based on the limited statistics at E > 2.5 TeV, the peak 
of the VHE emission related to HESS J 1023-575 appears to 
shift towards the location of the LAT pulsar PSRJ1022- 
5746 at higher energies although the statistic is too scarce 
to resolve further energy-dependent morphologies. The an- 
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gular separation of the two best-fitted positions for the 
low and high energy images is ^0.08°, of the order of the 
H.E.S.S. PSF. 



2.2. Spectral Analysis 

Spectral information has been obtained in the following 
way: Two circular regions were defined around the respec- 
tive centroids of HESS J1023-575 and HESS J1026-582, 
with radius 0.33°. The background was evaluated by us- 
ing the reflected background method, in which symmetric 
regions, not contamin ated by the sources, ar e used to ex- 
tract the background (jAharonian et al.ll2006bl ). The energy 
spectra are de rived by means of a forward-folding maximum 
likelihood fit (jPiron et al.ll2001l) . 

The two energy spectra are shown in Fig. [4] The spec- 
trum of the region containing HESS J1023-575 is well fit 
by a power law function dN/dE=N (E/l TeV)~ r , with a 
photon index of r=2.58±0.19 s t a t±0.2 sys and a normaliza- 
tion constant of N =(3.25±0.50 s tat)xl0- 12 TcV^cm-V 1 
(red markers and line), confirm ing the results from th e 
previous analysis (dashed line) (jAharonian et al.l l2007aT ). 
The previously measured flux level was slightly higher but 
this difference is explained by the larger integration radius 
(0.39°) used previously. At high energies the previous flux 
of HESS J1023-575 was contaminated with some fraction of 
the events related to the new source HESS Jf 026-582. The 
energy spectrum of HESS J1026-582 is also well represented 
by a power law function (blue markers and line), show- 



ing a hard photon index of r=1.94±0.20 sta t±0.2 sys as ex- 
pected from the energy maps, with a normalization con- 
stant of N =(0.99±0.34 stat )xl0- 12 TcV^cm^s- 1 . The 
systematic error on the normalization c onstant Np is es- 
timate d from simulated data to be 20% (jAharonian et al.l 
l2006bTl . 

The two measured photon indices in the energy range 
between 1 and 10 TeV differ by Ar=0.7, which taking into 
account their statistical errors, lead to an incompatibility 
of 2.5a for the two VHE emission regions. The different 
spectral indices as well as the clear separation of the two 
sources in the images favor the interpretation of the two 
emission regions as two independent sources rather than 
having a common origin. 



3. Discussion 

Follow up observation on the previously discovered source 
HESS J1023-575 lead to confirmation of the features of this 
source as well as the discovery of a new VHE 7-ray source, 
HESS J1026-582. The two TeV emission regions appear 
spatially distinct and show different spectral characteris- 
tics, favoring an interpretation of being two independent 
sources. Accordingly, different scenarios regarding the na- 
ture of both sources in the context of being related to ener- 
getic pulsars or massive stars and their winds are discussed 
in the following. 

3.1. Emission scenarios related to pulsars and their winds 

The spatial coincidence with recently discovered energetic 
7-ray pulsars motivates the investigation of scenarios where 
7-ray pulsars are powering a rclativistic wind nebula. In 
this scenario particles are accelerated to very high en- 



ergies along their propagation into the pulsar surround- 
ings or at the shocks produced in collisions of the winds 
with the surrounding medium. As a result of the interac- 
tions of rclativistic leptons with the magnetic field and low 
energy radiation (of synchrotron origin, thermal, or mi- 
crowave background), non-thermal radiation is produced 
from the lowest possible energies up to ^100 TeV. On the 
other hand, for magnetic fields of a few /iG, young elec- 
trons create a small synchrotron nebula around the pul- 
sar which should be visible in X-rays, in contrast to a 
much larger TeV nebula, generated by th e inverse Compton 
(IC) process (for a recent review see iGaensler Sz Slanel 
120061 ). Typically onl y pulsars with high spin-d own en- 
ergy (~10 33 erg s" 1 . [Kargaltsev fc Pavlov I l2010l ) produce 
prominent pulsar wind nebulae (PWNe). Both pulsars re- 
ported from Fermi LAT observations at energies below the 
H.E.S.S. detection threshold show high spin-down luminos- 
ity and features similar to other pulsars previously asso- 
ciated with VHE 7-ray sources such as HESS J1825-137 
(lAharonian et al.ll2006al) . 

3.1.1. HESS J1023-575 

The LAT-discovered 7-ray pulsar PSR J1022-5746 is the 
youngest and most energetic among the 7-ray only pul- 
sars known to date, being only one order of magnitude 
weaker than the Cr ab pulsar. An X-ray so urce, CXOU 
J102302.8-574606.9 (jTsuiimoto et al. I I2007D . found in a 
130 ksec Chandra observation has been suggested as a po- 
tential counterpart, 8' apart from the Westerlund 2 core, 
and its relation to the stellar cluster thus remains unknown. 
A distance to CXOU J102302.8-5 74606.9 of >10 kpc was 
propose d based on the column density from Chandra obser- 
vations (jDormodv et al.ll2009j) . The fitted column density 
Nh ^1.3x10 22 cm" 3 indicates a large distance d>8 kpc, 
although given the uncertainties in the measurement and 
the large absorption expected in the line of sight towards 
Westerlund 2, a possible association with the open cluster 
cannot be disregarded. On the other hand, no extended syn- 
chrotron PWN has been found yet albeit Suzaku observa- 
tion time was dedicate d to directly i nvestigate this hyp oth- 
esis (jFuiita et al.ll2009t) . As noted bv lFuiita et"aTl (|2009l) the 
ratio of TeV to X-ray luminosity is larger than one, assum- 
ing that the size of the X-ray emission they considered to 
obtain upper limits on the diffuse non-thermal emission in 
the 0.7 to 10 keV range corresponds to the same emitting 
region at VHE. That would suggest that the population 
of TeV emitting particle has already cooled and does not 
shine significantly in X ray syn chrotrons as is the case of 
most middle aged TeV nebulae (jde Jaeer fc Diannati-Ata'J 
120081) . 

With its high spin-down luminosity, the Fermi LAT pul- 
sar is sufficiently strong to power the observed VHE 7-ray 
emission. Assuming the pulsar is at a distance of 8 kpc, 
the luminosity of the H.E.S.S. source in the 1 to 10 TeV 
range is 4.8xl0 34 erg s _1 . The association with the pul- 
sar would imply conversion from rotational energy into 
non-thermal emission with efficiency 0.4%, comparable to 
those inferred for other VHE 7-r ay sources associated w ith 
PWNe, i.e., HES S J1912+1011 (lAlmronian et al.l l2008l) or 
HESS J1718-385 (jAharonian et al.l l2007bl ). However, the 
pulsar distance and age limit the maximum ex tension of 
the nebula (see i.e. iKennel fc Coronitil 1 1984T) and ren- 
der it difficult to explain the total extension of the TeV 
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Fig. 3: Profile of the TeV emission along the slice defined to 
contain the maximum of the two 7-ray excesses, as illustrated 
in the inset. The width of the slice is defined to be twice the 
H.E.S.S. angular resolution. The profile on the low energy image 
is shown in red while the one for the high energy image is shown 
in blue. The relative positions of WR 20a and the two Fermi 
LAT pulsars with respect to the center of the slice are indicated 
by dashed lines. 
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Fig. 4: Differential energy spectra of HESS J 1023-575 (red filled 
circles) and dash ed black line for the p revious detection and 
HESS J1026-582 (|Aharonian et alj|2007ej ) (blue open circles). 



emission as due to inverse Compton scattering of the 
PWN on the interstellar radiation fie ld (jGaensler &; Slanel 
120061 Ide Jaeer fc Diannati-Atall2008t ). The projected size 
of HESS J1023-575 corresponds to ~28 pc(d/8.0 kpc) (i. 
c. 6 or 31 pc assuming distances of 2 and 10 kpc, respec- 
tively). An estimation of the maximu m size of the nebula 
is given bv lKennel fc Coroniti I (|1984l) (KC model): 



R ~ 2 x 10 



2!! 
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^SD 



V 10 37 er£ 
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1/3 



( 1 



V10 5 yr 



1/3 



10 _12 erg/cm 2 



-1/3 



19 pc 



(1) 



obtained by balancing the total luminosity at present pro- 
vided by the pulsar against the interstellar medium pres- 
sure. The obtained value is 1.5 times smaller than the esti- 
mated size of the H.E.S.S. source for a pulsar as far as 8 kpc 
(28 pc) , although this value has a strong dependency on the 
assumed input parameters. Under these assumptions (using 



the pulsar parameters and a typical pressure for the ISM of 
p ism -2x10" 12 eV cm" 3 ), the KC PWN model for a pulsar 
located at the nominal distance of Westerlund 2 does not 
describe well the VHE emission and more sophisticated sce- 
nario are needed, including possible anisotropy of the pulsar 
wind or evolution of the spin-down luminosity which will 
allow a higher intrinsic spin-down energy. Also, the pres- 
ence of stellar winds from O type stars are known to create 
cavities, allowing a larger PWN size compared to expansion 
into more typical ISM pressures. A PWN according to the 
model of KC would be consistent with a distance estimate 
of ~2 kpc. Recently , a rev ised distance was estimated by 
ISaz Parkinson et all (|201Q[) suggesting PSRJ1022-5746 at 
a distance of 2.4 kpc. 

3.1.2. HESS J1026-582 

The association of the VHE 7-ray source HESSJ1026- 
582 with the Fermi LAT pulsar PSR J1028-5819 
(a=10 h 28 m 39.8 s and <5=-58°17'13", J2000, with a 
95% confidence level radius of 0.079°) seems likely and can 
be naturally described by interaction of the pulsar wind 
with the surrounding medium. Moreover, an alternative 
counterpart has not been proposed so far. The spin-down 
energy of the pulsar is high enough to power the VHE 7-ray 
emission, which, assuming a distance of 2.3 kpc, results 
in a total luminosity in the 1 to 10 TeV energy range of 
2.5 xlO 33 erg s _1 . The implied conversion efficiency from 
spin-down power into particle luminosity would imply a 
similar value to the one considered above, about 0.3%. 
No PWN association at X-ray energies has been found 
so far, although observations with both XMM-Newton 
and Suzaku have been carried out triggered by the Fermi 
LAT detection. The corresponding projected extension 
of HESS J1026-582 is 5.6 pc for an assumed distance 
of 2.3 kpc. Using Eq. [TJ the estimated maximum size is 
~8 pc, which is compatible with the observed extension at 
VHE energies. The large offset of the pulsar with respect 
to the center of gravity of the VHE source (0.28°) and 
asymmetric morphology is similar to other relic TeV PWN 
such as HESS ,11825-137 and can be explained either by the 
proper velocity of the pulsar or by an asymmetric reverse 
shock resulting from a supern ova that exploded initially 
in an inhomogeneous medium ([Blondin et al.ll2~001l) . If the 
pulsar was born at the H.E.S.S. source position, a pulsar 
transverse velocity of ^100 km s _1 can be estimated using 
the dispersion-measured based distance of 2.3 kpc, and a 
characteristic age of 8.9 xlO 4 yrs. This is in the observed 
range of known pulsar velocities. 

A summary of the observed parameters of the high en- 
ergy pulsars and inferred luminosity above 1 TeV for the 
nominal distances (d nom ) of the associated TeV sources is 
given in Table [T] 



3.2. Emission scenarios related to massive stars and their 
winds 

The possibility of an association of HESS J1023-575 with 
the open cluster Westerlund 2 was discussed in th e 
HESS J1023-575 discovery paper (|Aharonian et al.ll2007aD . 
Different mechanisms involving cosmic-rays accelerated in 
expanding stellar winds or supernova blast waves interact- 
ing with the boundaries of the interstellar radiation field or 
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Table 1: Properties of high energy pulsars associated to the TeV sources and inferred TeV luminosity for the nominal pulsar 
distance 



HESS source 


Pulsar Name 


E [erg s- 1 ] 


dnom [kpc] t [kyrs] 


B [10 12 G] 


L 7 (d/d nom ) 2 (l-10 TeV) [10 33 erg s" 1 ] 


HESS J1023-575 


PSRJ1022-5746 


l.lxlO 3 ' 


8 4.6 


6.6 


48 


HESS J1026-582 


PSRJ1028-5819 


8.3xl0 35 


2.3 89 


1.2 


2.5 



the boundaries of the radio blister ([Voelklll983|) are among 
the possible scenarios to explain the VHE 7-ray emission. 
The total mechanical luminosity (Wtotal) available from 
stellar wi nds in Westerlun d 2 is estimated to be 3.6 x 
10 51 erg (jRauw et al.ll2007t) . Obser vations of th e molecu- 
lar cloud content in the region by iDamd (|2007l ) shows a 
complex and rich environment and triggered high resolu- 
tion (1.5') observations with the Nanten 2 4-m sub-mm 
telescope. The latter observations resolved massive molec- 
ular clouds at ~16 km s" 1 and ^4 km s" 1 (6 . 5 and 5.2 
kpc), that were previously reported by IDamd (|2007f ). as 
well as a second component at ~-4 km s _1 which appears 
to be coincident with a region of ongoing star formation. To 
illustrate the extent of the star formation region, a three- 
color Spitzcr/IRAC GLIMPSE image is shown in Fig. [5j 
in red the 5.8 /xm image, 8 /zm in green and 3.6 /urn in 
blue. The Nanten 2 radio signal reported in the velocity 
ranges of 10.8 to 20.9 km s _1 and 1.2 to 8.7 km s^ 1 is over- 
laid in blue and purple dashed lines respectively, showing 
a strong spatial correlation with the cluster. The H.E.S.S. 
low energy significance contours are shown in white. These 
GMCs provide a sufficiently dense target (~8-9xl0 4 M Q 
which translates to a density of n^54-60 cm -3 for a size 
corresponding to the HESS J1023-575 extension at 8 kpc) 
for high-energy particle interactions, allowing the produc- 
tion of 7-rays from 7r-decay via inelastic pp-collisions. In 
both cases, considering the GMC mass and distance, (and 
using W PP =L 7 Tpp(n) , with t pp being the cooling time for 
p-p interactions) less than 0.05% of the total mechanical 
luminosity estimated for the stellar winds in Westerlund 2 
would be sufficient to explain the observed VHE emission. 

A power-law extrapolation of the H.E.S.S. spectrum to 
lower energies (to 30 GeV) yields a flux of 2.8x 10~ 8 TeV" 1 
enr 2 s _1 (or 1.7-6.2 xl0~ 8 TeV -1 cm~ 2 s _1 considering 
statistical errors in the normalization and spectral index). 
This flux extrapolation is well above the differential flux 
measured by Fermi LAT at ap proximately the s ame energy 
(~3.3xl0~ 9 TeV" 1 cm" 2 s'^ lAbdo et"alll2010h . Thus, the 
combination of the H.E.S.S. and Fermi LAT results im- 
plies a spectral hardening towards lower energies, which 
can be described with a turnover of the H.E.S.S. spectrum 
at ^100 GeV. If hadronic emission processes are respon- 
sible for the emission, the lower energy cutoff could result 
from cosmic-rays of energy <1 TeV being prevented from 
reaching t he GMC, perhaps via s low d iffusion as demon- 
strated by lAharonian fc Atovanl (fl996). where only high 
energy hadrons reach the surrounding molecular clouds. In 
this scenario, considering an age for Westerlund 2 of 1 to 
2 million years, low diffusion coefficients in the order of 
10 26 cm 2 s _1 are necessary to prevent low energy particles 
from reaching the 7-ray production site. 

Emission scenarios where the bulk of the VHE 7-rays 
is produced close to the massive WR stars WR20a and/or 
WR 20b were already disfavored given the observed exten- 
sion and lack of variability in t he VH E 7-ray emission, 
as discussed in lAharonian et all (l2007a|h However, 7-rays 




Galactic Longitude (deg) 



Fig. 5: Spitzer/IRAC GLIMPSE three-color image (in red the 
5.8 /jm image, 8 /im in green and 3.6/xm in blue) of RCW 49 
overlaid with 12 CO(J = 2-1) contours at velocity ranges of 10.8 
to 20.9 km s _1 and 1.2 to 8.7 km s _1 observed with the Nanten 2 
telescope (in blue and purple respectively). The white contours 
show the H.E.S.S. significance derived for a correlation radius of 
0.1°, hard cuts and a maximum energy of 2.5 TeV. The lowest 
contour corresponds to 4 a. 



from collective stellar wind interactions, where energetic 
particles experience multiple shocks, may still contribute 
to the observed VHE emission. 

4. Conclusions 

The H.E.S.S. source HESS J1023-575 has been re-observed 
in order to search for energy-dependent morphology, that 
could provide clues on the origin of the 7-ray emission. The 
analysis of the new data unveils a second VHE emission 
peak towards the direction of the energetic high energy 
pulsar PSR J1028-5819, reported recently from Fermi LAT 
observations at lower energies (0.1 to ~30 GeV). The cen- 
ters of gravity of the two VHE 7-ray emission regions are 
well separated (0.6°) and the radial profile shows two clear 
peaks at the higher energies, while at lower energies the 
VHE emission corresponding to HESS J1026-582 appears 
very much suppressed. This new source is detected with a 
significance level of 7tr and is characterized by a very hard 
photon spectral index of 1.9. HESS J1026-582 is preferably 
interpreted as due to a PWN energized by PSR J1028-5819, 
with similar characteristics as those previously detected by 
H.E.S.S. (such as hard spectrum photon index, offset of the 
high energy pulsar with respect the center of the TeV emis- 
sion, similar conversion efficiency from rotational energy 
into non-thermal emission, etc). Further observations with 
X-ray satellites of this region will have a chance to provide 
the required additional information to support and settle 
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the suggested scenario of being another GeV Pulsar/TeV 
PWN association. 

The previously reported VHE source HESS J1023-575 
has been confirmed and the analysis of the new data shows 
compatible results. Multi- wavelength observations of the re- 
gions provided new insights on the potential origin of the 7- 
ray emission as discussed in the previous section. A possible 
association with the recently discovered 7-ray PSRJ1022- 
5746 (1FGL J1023. 0-5746), is supported by spatial coinci- 
dence with HESS J1023-575, and can be easily accommo- 
dated in terms of energetics if the pulsar is as close as 2-3 
kpc, while more sophisticated PWN models are necessary if 
it is located significantly further away. On the other hand, 
the molecular content of the region provides sufficient tar- 
get material to explain the emission through hadronic in- 
teraction of cosmic-rays accelerated by a range of feasible 
mechanisms in the open cluster interacting with molecular 
clouds. 

Whether the emission is due to the not-yet-detected 
PWN associated to the Fermi LAT pulsar PSR J1022-5746 
or mechanisms related to acceleration of cosmic-rays in the 
open cluster, a spectral flattening at lower energies must 
be present to explain the lower MeV-GeV flux level. In 
the first case, IC processes would naturally reproduce the 
spectral shape, while for hadronic emission, a lower cutoff 
at a several tens of GeV would reflect a lower cutoff on 
the cosmic-ray spectrum at ~ 1 TeV. The latter scenario 
could be explained by assuming a low diffusion coefficient 
in the region, which would not be surprising in the context 
of the cavities created by the stellar winds from the young 
OB stars in the stellar cluster. Despite the already accom- 
plished multi-wavelength coverage, a firm identification is 
still pending, whereas the high-power pulsar PSRJ1022- 
5746 is the prime counterpart hypothesis to adopt now. 
Unfortunately the statistics at VHE is too scarce to dis- 
tinguish any further hint of energy- dependent morphology 
which would have the potential to settle the association 
with the pulsar or the open cluster. A clear identification 
of a PWN associated to the Fermi LAT pulsar PSR J1022- 
5746 in X-rays would help to elucidate the real nature of 
HESS J1023-575, but given the stellar density and richness 
in X-ray sources towards Westerlund 2, this appears to be a 
serious observational challenge. In particular deeper VHE 
observations might allow higher energy resolution studies 
and a firm identification of the origin of the VHE emission. 
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